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Multi Solitons of a Bose-Einstein Condensate 
in a Three-Dimensional Ring 
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A Bose-Einstein Condensate (BEC) made of alkali-metal atoms at ultra-low temperatures is well de- 
scribed by the three-dimensional cubic nonlinear Schrodinger equation, the so-called Gross-Pitaevskii equa- 
tion (GPE). Here we consider an attractive BEC in a ring and by solving the GPE we predict the existence of 
bright solitons with single and multi-peaks, showing that they have a limited domain of dynamical stability. 
We also discuss finite-temperature effects on the transition from the uniform phase to the localized phase. 
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Single and multiple bright soliton configurations 
in a Bose Einstein condensates (BEC) have been 
observed in two different experiments ^ |2] ■ These 
solitons, which can travel for long distances with- 
out dispersion, have been theoretically investigated 
by various authors due to their relevance in non- 
linear atom optics [Hll^- Very interesting appears 
the not yet experimentally studied case of an at- 
tractive BEC in a ring. For this system it has been 
predicted, on the basis of the one-dimensional (ID) 
theory, a quantum phase transition from a uniform 
condensate to a bright soliton j5j . 

Here we consider an attractive BEC in a 3D ring, 
taking into account transverse variations of the BEC 
width. We show that the phase diagram of the sys- 
tem reveals novel and peculiar phenomena: the lo- 
calized soliton has a limited existence and stabil- 
ity domain and the system supports also multi-peak 
solitons, which can be dynamically stable. 

The time-dependent Gross-Pitaevskii equation 
(GPE) for a BEC in a trapping potential U{t) is 
given by 



where is the number of condensed atoms and 
as the s-wave scattering length of inter-atomic in- 
teraction. The GPE describes A^ trapped bosons 
in the same quantum state ^(r, i): a pure BEC 
of A^ atoms. The GPE is a mean- field Hartree 
equation and it is accurate for a dilute gas of con- 
densed bosons at ultra-low temperature. The cubic 
term is obtained by using a contact pseudo-potential 
V{y,y') = Airh'^as/m 6^{r — r') to model the two- 
body interaction. If > there is repulsion be- 
tween particles, if < there is attraction between 
particles [S]. 

We describe an attractive BEC in a cylinder, i.e. 
under transverse radial harmonic confinement and 
without axial longitudinal confinement, by using 
the self-focusing GPE (og < 0) with u;± the trans- 
verse frequency of harmonic confinement. This self- 
focusing 3D GPE is the Euler-Lagrange equation of 
the following Lagrangian density 
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where T = 4:Tt{N — l)\as\/a± is the scaled inter- 
atomic strength. Here we use scaled units: energy 
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in units ?ilo±, time in units of wj^, length in units 
of a± = ^/W/JnuJ±)- To simplify the 3D prob- 
lem, let us choose the trial wave function ■0(r, t) = 
/(z,t)exp [-(x^ +y2)/2o-(z,t)2]/vrcr(z,t)2. Insert- 
ing it into the Lagrangian, integrating over x and y 
and neglecting the derivatives of a{z, t) one finds an 
effective Lagrangian density C which depends on the 
axial wave function f{z, t) and the transverse width 
a{z,t). The two Euler-Lagrange equations of the 
effective Lagrangian C with respect to / and a are 



f 



/ 



l-9\f\' 



where g = T/{2tt) = 2{N - l)\as\/a_L. For a = 1 
the first equation becomes the familiar self-focusing 
ID GPE, also called cubic nonlinear Schrodinger 
equation. The equation for a can be inserted 
into the equation for /. In this way one gets the 
time-dependent non-polynomial Schrodinger equa- 
tion (NPSE) 0. 

In 7 we have verified that the NPSE gives re- 
sults always very close to that of the 3D GPE un- 
der transverse harmonic confinement for both pos- 
itive and negative scattering length a^. The self- 
focusing NPSE (os < 0) admits a shape-invariant 
solution, called bright soliton, due to the inter- 
play between the dispersive kinetic term and the 
attractive interaction term. Note that bright soli- 
tons of the self-focusing 3D GPE do not exist with- 
out transverse confinement. We have found ^ 
that the bright soliton reduces, in the limit of very 
strong transverse confinement, namely for g « 1 
(ID regime), to the text-book bright soliton of the 
self-focusing ID GPE, whose axial wave function 
is f{z) = \fgT^ sech{gz/2) Contrary to the self- 
focusing ID GPE, the bright soliton of the self- 
focusing NPSE has a limited domain of existence: 
< 5 < |. If the strength g exceeds 4/3 there is the 
so-called collapse of the BEC. Thus, the transverse 
dynamics of the BEC induces the instability of the 
solitonic solution. 

As previously stressed, bright solitons have been 
observed in two experiments: a single soliton in ex- 
pulsive trap at the Ecole Normale Superieure (ENS) 



of Paris |2I, and a soliton train in a harmonic trap 
at Rice University jT]. We belive that in the next 
future there will be experiments on attractive BEC 
in a toroidal trap. In fact, toroidal traps have been 
recently produced: at Georgia Tech cold atoms in a 
magnetic ring (diameter of 2 cm) [5] and at Berke- 
ley a BEC in a magnetic ring (diameter of 3 mm) 

To model an attractive BEC in a toroidal trap 
we consider now the self-focusing NPSE where the 
solution f{z,t) has periodic boundary conditions: 

f{z + L,t) = f{z,t) 

where L = 2ttR is the azimuthal longitudinal length 
of the BEC and R is the azimuthal radius of the 
trap. Due to the periodic boundary conditions the 
stationary self-focusing NPSE admits a uniform so- 
lution, given by f{z) = Xj^fL. One sees immedi- 
ately that the uniform solution exists only for g < 
2ttR. The stationary NPSE admits also nodeless 
and nodal localized solutions (solitonic solutions). 
In Fig. 1 we plot the existence diagram of these 
solitonic solutions in the {R, g) plane. The diagram 
is obtained by numerically solving the NPSE. The 
one-peak localized nodeless solution exists between 
the two solid lines of Fig. 1. Nodeless two-peak lo- 
calized solution exists between the two dashed lines 
while the nodal two-peak localized solution exists 
between the two dot-dashed lines. 

We analyze the dynamical stability of the so- 
lutions by numerically solving the time-dependent 
NPSE taking as initial condition the stationary lo- 
calized solution f{z) with a very weak perturbation. 
In the left panels of Fig. 2 we plot the density profile 
p{z) = |/(2)p of the solitonic solutions. In the right 
panels we show the time-dependence of the mean 
squared widths (z^) for the weakly perturbed soli- 
tonic solutions. A periodic oscillation means that 
the solitonic solution is dynamically stable. In the 
calculations we choose a ring axial lenght L = 15 
and an interaction strength g = 1. We have stud- 
ied the dynamical stability for various initial con- 
ditions. The main results are the following: i) the 
uniform solution is stable only below the lower curve 
of existence of the one-peak soliton; ii) the one-peak 
soliton is stable where it exists; iii) the nodal two- 
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FIG. 1. Existence diagram in the (i?, g) plane, where g — 
2A^|as|/aj_ is the interaction coupling and B = L/ (27r) is 
the azimuthal radius of the ring (in units a±). Uniform 
solution: exists for g < 2ttR. One-peak localized node- 
less solution: exists between the two solid lines. Node- 
less two-peak localized solution exists between the two 
dashed lines. The nodal two-peak localized solution ex- 
ists between the two dot-dashed lines. 
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FIG. 2. Left panels: density profile p{z) of the soli- 
tonic solutions. Right panels: time-dependence of the 
mean squared widths {z^) for the weakly perturbed soli- 
tonic solutions. Ring axial lenght: L = 15. Interaction 
strength: g = l- 
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peak soliton is stable in the plane (i?, g) only below 
the upper curve of existence of the one-peak soliton; 
iv) the nodeless two-peak soliton is unstable. Sim- 
ilar results are found for solitonic solutions with a 
larger number of peaks. 

To conclude, we observe that the critical strength 
g* = TTa±/2R of the transition from the uniform 
state to the localized state at zero temperature has 
been determined in [5] by investigating the Bogoli- 
ubov energy spectrum of the the self-focusing ID 
GPE. The inclusion of finite-temperature effects is 
quite simple if one consider a non-interacting ther- 
mal cloud. In this way we find that the critical 
temperature T* of transition from the uniform to 
the localized state is given by 

huJl, 4i?2 4i?|a^| j 

where /c^ is the Boltzmann constant, and using not- 
scaled quantities. Obviously if < TTa\/ {AR\as\) 
the system remains uniform at all temperatures. 
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